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ABSTRACT 

The discovery of an optical counterpart to GRB0 10222 (detected by BeppoSAX; Piro 2001) was announced 4.4 
hrs after the burst by Henden (2001a). The Sloan Digital Sky Survey's 0.5m photometric telescope (PT) and 2.5m 
survey telescope were used to observe the afterglow of GRB010222 starting 4.8 hours after the GRB. The 0.5m 
PT observed the afterglow in five, 300 sec g* band exposures over the course of half an hour, measuring a temporal 
decay rate in this short period of F v oc ? _I 0±0 - 5 . The 2.5m camera imaged the counterpart nearly simultaneously 
in five filters (u*g*r*i*z*), with r* = 18.74 ±0.02 at 12:10 UT. These multicolor observations, corrected for 
reddening and the afterglow's temporal decay, are well fit by the power-law F v oc v~ {) ' 90±0 03 with the exception of 
the u* band UV flux which is 20% below this slope. We examine possible interpretations of this spectral shape, 
including source extinction in a star forming region. 

Subject headings: gamma rays: bursts (GRB010222) 



1. introduction 

Gamma Ray Bursts (GRBs) were first detected over three 
decades ago by the Vela satellites (Klebesadel et al. 1973), 
and the first search for optical counterparts started nearly im- 
mediately with W. A. Wheaton's use of the Prairie Network 
(Grindlay & Wright 1974). These searches were fruitless un- 
til very recently; positions accurate to a few arcminutes were 

' not available for days, after the bursts had decayed substan- 
tially, placing afterglows beyond the reach of the few large tele- 
scopes searching for them. BATSE's near-real-time coordinates 
had several degree positional errors (Paciesas et al. 1999) al- 
lowing only specialized wide field instruments to respond to 

■ its triggers (Krimm et al. 1995; Lee et al. 1997; Akerlof et al. 
1999). The BeppoSAX satellite (Scarsi 1993) was the first to 
provide arc minute accuracy within a few hours of a GRB. With 
the early announcements of those accurate positions, beginning 
in 1997, large telescopes could join the search and discovered 
GRB afterglows starting with GRB970228 (Groot et al. 1997; 
van Paradijs et al. 1997). 

The following work describes observations of GRB0 10222 
with the Sloan Digital Sky Survey's (SDSS; York et al. 2000) 



telescopes. The SDSS is a project to image 10,000 deg 2 of the 
Northern Galactic Cap in five different filters (u*,g*,r*,i*,z*) 
to a depth of r* ~ 23 and to perform followup spectroscopy of 
the 10 6 brightest galaxies and 10 quasars found in the photom- 
etry. The SDSS is designed to be on the u 1 g'r'i'z 1 photomet- 
ric system described in Fukugita et al. (1996) which is an AB V 
system where flat spectrum objects (F v oc v°) have zero colors 
(Fukugita et al. 1996). The magnitudes in this paper are quoted 
on the preliminary u* ,g* ,r* ,i* ,z* system which may differ by 
at most a few percent from the system of Fukugita et al. (1996). 
The dedicated survey instruments, a 2.5m survey telescope and 
a 0.5m photometric telescope (PT), are located at Apache Point 
Observatory (APO) in Sunspot, New Mexico. 

2. observations 



GRB0 10222 was detected by BeppoSAX on 2001 Febru- 
ary 22 at 07:23:30 U.T in both the Gamma Ray Burst Mon- 
itor (GRBM; 40-700 keV) and the Wide Field Camera Unit 1 
(WFC1; 2-28 keV) instruments, and was "possibly the brightest 
(GRB) ever observed by BeppoSAX" (Piro 2001). The coordi- 
nates of the BeppoSAX detection were distributed via the GRB 
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Coordinates Network (GCN; Barthelmy et al. 1998) at 10:36 
UT (Piro 2001), and Henden (2001a) reported the discovery of 
an optical counterpart at 1 1 :48 UT, 4.4 hours after the trigger. 
(See finding chart, Figure 1.) At this time conditions at APO 
were not ideal for SDSS survey imaging as clouds were ap- 
proaching, and the time remaining in the night did not allow 
for a switch to spectroscopy; thus SDSS observers decided to 
follow up the counterpart with both the 0.5m PT and the 2.5m 
survey telescope. Fortunately for these observations the cloud 
passed before GRB imaging began and conditions were more 
photometric after the cloud than before. 

2.1. 0.5m Photometric Telescope Observations 

The photometric telescope is an f/8.8, 0.5m telescope 
equipped with u*g*r*i*z* filters. The single SITe 2048 x 2048 
CCD camera has a 41.5' x 41.5' field of view. The PT took 
a series of five, 300 second observations in g* band centered 
on the reported GRB0 10222 location, following the afterglow 
for approximately 30 minutes before ending operations for 
the night (see Table 1). Normally, the photometric telescope 
and the associated reduction software are used on objects with 
g* < 18.0. Since the GRB exposures were unusually long and 
the counterpart was relatively dim (the Poisson error limit is 
~ 3%) photometry was performed within a smaller than stan- 
dard aperture to improve the relative photometry of faint ob- 
jects. The counterpart magnitudes were then corrected using 
a sigma clipped mean of the magnitude offsets in each frame 
from the mean magnitude across the five frames of well mea- 
sured stars (g* < 17.0); corrections were at most 0.005 mag, in- 
dicating conditions were photometric over the 30 minute times- 
pan. 

2.2. 2.5m Survey Telescope Observations 

The 2.5m survey telescope is an f/5, 3° field of view tele- 
scope designed and constructed for the SDSS. The telescope 
has two interchangeable instruments, an imaging camera and a 
fiber-fed spectrograph. The imaging camera (Gunn et al. 1998) 
includes an array of 30 2048 x 2048 CCDs in six columns of 
five CCDs each, one CCD for each of the 5 filters. The camera 
operates in a drift scan mode, scanning the sky in great circles 
at sidereal rate. Astronomical objects are imaged for 53.9 sec- 
onds in each CCD in the order r*i*u*z*g*. Because of the gaps 
between columns the telescope must observe a second such in- 
terleaved strip to make a complete stripe. 

For GRB010222 the 2.5m telescope observed two short in- 
terleaved strips covering a roughly 2.5° square region. The 
GRB0 10222 counterpart was found in the second strip, field 
22 of camera column 3, run 2143. The images were pro- 
cessed through the normal SDSS data processing pipelines and 
calibrated against two 0.5m PT secondary fields (hereafter, 
patches) of the GRB field observed March 14 and 17. These 
patches in turn were calibrated against a system of standard 
stars (Smith et al., in preparation) which the 0.5m PT observes 
several times throughout the night to measure extinction and de- 
termine photometric zeropoints. Our diagnostic tests of the lo- 
cation of the stellar locus and number counts of various classes 
of objects, as compared with the approximately 1000 deg 2 of 
sky observed already in the survey, indicate the relative calibra- 
tions are no worse than 2% in any filter. The zeropoints also 
agree to within 1 % with a second indirect calibration based on 
other 2.5m data from the same night and four secondary patches 
from two previous nights. Thus we are confident the relative 



(absolute) errors are no greater than the standard SDSS val- 
ues of 3%(5%) for u*, g*, and z* , and 2%(3%) for r* and i* 
(where the absolute errors include possible differences between 
the SDSS preliminary photometric system and the system of 
Fukugita et al. 1996). The calibrated 2.5m magnitudes are 
shown in Table 1 along with 0.5m PT observations. Table 2 in- 
cludes u*g*r*i*z* magnitudes for eleven reference stars in field 
22, including the reference star "A" listed in McDowell et al. 
(2001). All stars have been selected from unsaturated and non- 
interpolated (for cosmic ray correction) stars with r* < 19.0 
with the exception of the reference star "A" from McDowell et 
al. (2001). This star was interpolated in r* band, but the correc- 
tion was within the quoted errors and visual inspection revealed 
no problems. 

3. RESULTS 

The fading of GRB optical afterglows is often well fit by a 
power-law decay of F v oc t a with a m -1 although decay rates 
from slightly less than this to a = -2.1 (Groot et al. 1998) and 
breaks to steeper power-laws have been observed in some after- 
glows. GRB010222 appears to be best fit by broken power-law 
models (Holland et al. 2001; Masetti et al. 2001) with early de- 
cay rates of a m -0.6 to -0.7 and steeper later time decays of 
a w — 1.3 with the break occurring around 0.5 days after the 
GRB. Considering the short time span and the limiting errors, 
measurements of the decay rate with the 0.5m PT observations 
were difficult. The error weighted least squares fit to the five g* 
points gives a = -1.0 ±0.5 (see Figure 2), which is consistent 
with other early decay rates measured for this burst (Price et 
al. 2001; Fynbo et al. 2001; Stanek et al. 2001b; Holland et al. 
2001; Masetti et al. 2001). 

The spectral shape of GRB afterglows is also well fit by a 
power-law, F v oc v 13 with typical values of « — 1. In or- 
der to derive a power-law fit for our 2.5m u*g*r*i*z* observa- 
tions we first corrected for the local Galactic extinction with 
the dust map of Schlegel, Finkbeiner, & Davis (1998) which 
gives extinction values of A u * = 0.118, A g * = 0.087, A r * = 
0.063, A,» = 0.048, and A z » = 0.034 at the reported location 
a = 14 h 52 m 12: s 55, S = +43°01'06"2 (J2000; Henden & Vrba 
2001b). We next wished to correct for the small effect of fading 
over the few minutes between exposures in the individual bands 
at 4.8 hours after the burst. Because of the large errors in our 
own decay rate measurement we instead used a least squares fit 
to a single power-law for all reported R band data points within 
8 hours after the burst (and before the w 0.5 day break; see Mc- 
Dowell et al. 2001; Stanek et al. 2001a; Watanabe et al. 2001; 
Holland et al. 2001) with magnitudes adjusted to the calibra- 
tion of Henden (2001c) where necessary. The resulting fit to 
the published R band data is a = -0.71 ± 0. 10, consistent with 
values reported by Holland et al. (2001). We applied decay cor- 
rections relative to r* of [-0.0035, -0.0066, -0.010, -0.013] to 
i*,u*,z*,g*. 

Once these corrections were applied, the least squares fit to 
all five filters is = -1.10 ±0.10. However, a much better fit 
can be obtained by excluding the u* filter; the remaining non- 
UV filters have a least squares fit of = -0.90 ± 0.03 (see Fig- 
ure 3). This second value agrees closely with the spectral slope 
of = -0.89 ± 0.03 observed by Jha et al. (2001b) in a spec- 
trum observed 4.92 hours after the burst, shortly after the 2.5m 
observations. Our u* magnitude, at an effective wavelength of 
3565A, is approximately 20% lower than the power-law fit to 
the other bands. A similar deficit was seen in the U band ob- 
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servations of Masetti et al. (2001) one and two days after the 
burst, and Jha et al. (2001b) may see the beginning of this steep- 
ening in their last binned spectra point near 4000A. (After this 
work was submitted similar U band results were reported by 
Stanek et al. (2001c) for observations from 0.2 to 2 days after 
the burst.) These independent observations indicate this spec- 
tral feature remained constant for at least two days. 

We propose that the break in the spectrum at u* may be an 
indication of one of two possibilities: either the Lya forest or 
extinction at the source. The first possible explanation for the 
u* deficit is that the counterpart is at a redshift near 2 rather than 
at the redshift 1 .477 absorption system reported by Garnavich 
et al. (2001); Jha et al. (2001a); Bloom et al. (2001); Castro et 
al. (2001); Jha et al. (2001b) and Masetti et al. (2001). For an 
object without detected emission lines such as GRB010222, an 
absorption line system can only provide a lower z limit to the 
source redshift. The spectrum of Jha et al. (2001b) does not ex- 
tend far past 4000A, thus no upper limit is imposed until z = 2.3. 
For observed QSOs the Lya forest enters u* at redshifts slightly 
above z = 1 .477, but the 20% depression observed here would 
not occur unless the redshift were > 2.0 (Cristiani et al. 1993). 
However Jha et al. (2001b) convincingly argue that given the 
strength of the z = 1 .477 absorber the GRB source is almost 
certainly at that redshift. Thus a true GRB source redshift of 
« 2 would seem to be an unlikely explanation. 

A second, more probable, explanation is that the counterpart 
may reside in a star-forming region at z = 1 .477 similar to the 
Large Magellanic Cloud (LMC) or Small Magellanic Cloud 
(SMC) and be extincted at the source. Dust in front of the 
GRB could cause the extinction of the afterglow and the gas 
would explain the large equivalent widths in the z = 1 .477 ab- 
sorption system (York et al. 1986). To examine this possibility 
we have fit the full extinction curve model of Reichart (2001) 
to the SDSS data. Acceptable fits can be found for a wide range 
of intrinsic power-law spectra (with index (3'). We present two 
possibilities, (3' = -0.75 and -0.5. These choices for f3' are mo- 
tivated as follows: under the models of Sari, Piran, & Narayan 
(1998) and Sari, Piran, & Halpern (1999) the afterglow is de- 
scribed in terms of synchrotron emission from a decelerating 
relativistic shell or jet colliding with the surrounding ISM. The 
resulting spectrum can be expressed as four power-laws bro- 
ken at three time-dependent frequencies, the synchrotron self- 
absorption frequency v a , the cooling frequency v c , and the fre- 
quency corresponding to the minimum Lorentz factor of accel- 
erated electrons v m . If the shock evolves adiabatically in a con- 
stant density medium, the break in the light curve at « 0.5 days 
(Holland et al. 2001 ; Masetti et al. 2001) might be explained by 
a jet if the observed v opt > v c and v m , and (3' s=a -0.75. If the 
shock instead evolves radiatively, the break in the light curve 
might be explained by v m passing through the optical at « 0.5 
days if v c < v opt < u m , and f3' « -0.5. For (3 1 = -0.75(-0.5), 
we find that the best fits in the Reichart (2001) model are the 
source extinction Ay ~ 0.032 (0.13) mag, the slope of the UV 
linear component C2 ~ 1.35 (1.34), the strength of the UV bump 
C3 « 8.1 (2.7), and the strength of the FUV non-linear compo- 
nent C4 « 30 (6.9). These curves are shown in Figure 4 which 
includes a typical SMC-like extinction curve for a source spec- 
trum with (3' = -0.75 and source extinction Ay = 0.10. For the 
f3' = -0.75 case there is a strong degeneracy between Ay and 
the parameters C3 and q such that only Ay ■ C3 and Ay ■ C4 can 
be constrained; C3 and C4 can be increased to any value by de- 
creasing Ay, thus statistically we can only set lower bounds. For 



(3' = -0.75, A v < 0.057 mag (la) and Ay > at the 4.8er confi- 
dence level; C2 = 1.35^;|f ; C3 > at the 1.1 a confidence level; 
and C4 > ll(lcr), C4 > at the 2.5a confidence level. Further, 
C4 > 1 (higher than any observed value) at the 2.1<r confidence 
level. For f3' = -0.5 the degeneracy is not as stong, and we can 
place the following limits: Ay = 0.1 3±g {j® mag and Ay > at the 
4.8cr confidence level; c 2 = 1-34+J];^; c 3 > at the 1.3er confi- 
dence level and C3 < 5.6 (la); c\ > at the 2.1cr confidence 
level, C4 > 1 at the 1.6cr confidence level, and c\ > 3.4 (la). 
The Reichart (2001) fits to the two GRB models are approxi- 
mately equally likely (with the /?' = -0.75 model fit only 1.4 
times more probable than the (3' = -0.5 model fit). 

The Reichart (2001) best fit values of C2 and the second value 
of C3 are typical of that observed in the LMC. Given the errors, 
the first value of C3 is not inconsistent with this interpretation. 
However, for both afterglow models the values of q, required 
to extinguish u* relative to the other bands, are about an order 
of magnitude greater than those found in the LMC or SMC. 
Waxman & Draine (2000) and Galama & Wijers (2001) pro- 
pose that the optical flash (e.g., Akerlof et al. 1999) and the 
burst may sublimate and fragment dust in the circumburst envi- 
ronment. If small (radius < 300 ) graphite grains, which may 
be responsible for the FUV non-linear component (Draine & 
Lee 1984), survive in greater numbers in this environment, this 
value of C4 is not unreasonable. Alternatively, the large C4 value 
and u* band deficit could be due to absorption by molecular 
hydrogen (Draine 2000) which would span the entire u* band 
at z = 1.477; however the expected feature at A < 1650A red- 
shifted to A m 4000A is not obvious in the published spectra. 

4. CONCLUSIONS 

The serendipitous 2.5m survey telescope observations of 
GRB0 10222 occurred in this case because of the very fortunate 
timing of a counterpart discovery announcement towards the 
end of a night when conditions did not favor normal survey op- 
erations. The 2.5m camera is an unwieldy instrument for rapid 
followup observations; nonetheless this observation has shown 
the value of early five filter observations. In addition to the 
measurement of the g*r*i*z* spectral slope (F v oc jy-°- 90 ±0 03^ 
the break to a steeper slope in u* (also seen in the U band ob- 
servations of Masetti et al. 2001 and Stanek et al. 2001c) was 
not predicted or seen in spectra, and may indicate an alternate 
source redshift, source extinction in a star forming region mod- 
ified by the GRB or its progenitor, or something else entirely. 

The 0.5m PT is an automated telescope and in general much 
better suited for GRB followup observations than the 2.5m sur- 
vey telescope. In this case the same timing that was so fortunate 
for the 2.5m was disadvantageous for the PT, which was only 
able to observe the burst near its limit and for a short period 
before shutting down for the night. Due to an afterglow's rapid 
decay, typical BeppoSAX delays of several hours place after- 
glows near the detection limit of smaller telescopes. HETE-2, 
launched in October of 2000, will soon provide ^10 arcminute 
positions for GRBs within minutes of the trigger, potentially al- 
lowing telescopes such as the PT to measure both the spectral 
and temporal behavior of a burst in the first few hours. 
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Table 1 



SDSS Observations of the GRB010222 Afterglow. 



UTC a 


telescope 


band 


exposure (sec) 


magnitude b 


12:09:35 


2.5m 


r* 


54 


18.74 ±0.02 


12:10:47 


2.5m 


i* 


54 


18.53 ±0.02 


12:11:59 


2.5m 


u* 


54 


19.56 ±0.03 


12:13:10 


2.5m 


z* 


54 


18.34±0.03 


12:14:22 


2.5m 


8* 


54 


19.02 ±0.02 


12:13:15 


0.5m 


8* 


300 


19.04 ±0.04 


12:19:36 


0.5m 


8* 


300 


19.05 ±0.04 


12:25:58 


0.5m 


8* 


300 


19.03 ±0.04 


12:32:19 


0.5m 


8* 


300 


19.10±0.04 


12:38:41 


0.5m 


8* 


300 


19.14±0.04 



a Exposure start time, 2001 Feb 22. 

b Statistical errors; absolute photometry errors for the 2.5m may 
be as large as 5% for u*, g*, and z*, 3% for r* and i*. 
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a» 


S b 




g* 


r* 


i* 


z* 


grb 


14 h 52 m 12*51 


+43°01'06"2 


19.56 ±0.03 


19.02 ±0.02 


18.74 ±0.02 


18.53 ±0.02 


18.34 ±0.03 


A 


14 h 52 m 07 s 51 


+42°58'48"6 


19.39 ±0.02 


17.96 ±0.02 


17.40 ±0.02 


17.18 ±0.01 


17.08 ±0.01 


B 


14 h 52 m 12*57 


+42°55'59"3 


18.29 ±0.02 


17.02 ± 0.02 


16.48 ±0.02 


16.26 ±0.01 


16.10±0.01 


c 


14 h 52 m 16 s 06 


+43°02'38"5 


20.02 ±0.03 


19.03 ±0.02 


18.72 ±0.02 


18.61 ±0.01 


18.56 ±0.03 


D 


14 h 52 m 21?86 


+42°56'29"0 


19.78 ±0.03 


18.42 ±0.02 


17.87 ±0.02 


17.67 ±0.01 


17.53 ±0.02 


E 


14 h 52 m 28 s 65 


+43° 02' 32" 6 


20.37 ±0.04 


18.50 ±0.02 


17.75 ±0.02 


17.45 ±0.01 


17.29 ±0.02 


f 


14 h 52 m 31?09 


+43°03'14"2 


17.06 ± 0.01 


15.89 ±0.02 


15.48 ±0.01 


15.30±0.01 


15.26 ±0.01 


G 


14 h 52 m 41?98 


+43°03'16"8 


19.37 ±0.02 


18.15 ±0.02 


17.69 ±0.02 


17.52 ±0.01 


17.42 ± 0.02 


H 


14 h 52 m 45 s 90 


+42°57'09"2 


20.07 ± 0.04 


17.80 ±0.02 


16.81 ±0.02 


16.46 ±0.01 


16.24 ±0.01 




14 h 52 m 50 s 13 


+42°55'22"7 


17.84 ±0.01 


16.51 ±0.02 


15.99 ±0.01 


15.83 ±0.01 


15.73 ±0.01 


J 


14 h 52 m 51 s 32 


+42°54'56"3 


18.53 ±0.02 


16.49 ±0.02 


15.63 ±0.01 


15.37 ±0.01 


15.18 ±0.01 


K 


14 h 52 m 5L85 


+43°03'05"7 


18.59 ±0.02 


16.91 ±0.02 


16.28 ±0.01 


16.01 ±0.01 


15.89 ±0.01 



^Selected from stars in field 22 with magnitude r* < 19.0 which were well measured, unsaturated, and non-interpolated (cosmic 
ray corrected) in all five filters, with the exception of the first star (A) in the list which is the reference star 'A" listed in McDowell et 
al. (2001). This star was interpolated in r* band but manual inspection revealed no problems; the correction was small and within the 
quoted errors. Listed errors are statistical only, absolute photometry errors for the 2.5m may be as large as 5% for «* , g* , and z* , 3% 
for r* and i* . 

b Because of the non-standard orientation and short length of this stripe the astrometric errors are unusually large, approximately 
0"3. 
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FIG. 1 . — GRB010222 2.5m telescope r* image. The image is 8 arcmin square with 1 arcmin tick marks. North is approximately 3° clockwise from up, east is to 
the left. Stars from Table 1 within this subsection of field 22 are indicated, including the reference star "A" of McDowell et al. (2001). 
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FIG. 2. — Relative photometry for five 0.5m PT g* band observations, along with the single 2.5m g* band observation. The best fit decay curve of the form 
F v oc t a to the five PT g* band observations is a = -1 .0 ± 0.5. 
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1-9 
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O 
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2.2 - 

2.1 - 
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FIG. 3. — The 2.5m multiband observations at a single epoch. The best fit to F v oc u@ with all five bands is j3 = — 1.10 ± 0.10, shown as a dotted line. Excluding 
u* band produces a fit of fi = -0.90 ± 0.03, shown above as a dashed line. 
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FIG. 4. — Extinction model fits to the 2.5m multiband observations. The two best fit extinction models from Reichart (2001) for fi' = —0.75 and 0.5 as discussed 
in section 3, as well as an SMC curve for an unextinguished spectrum with /3' = —0.75, are presented here. The green and blue curves through the data points 
correspond to the Reichart (2001) extinction model fit for z = I ATI and a source spectrum of fi' = —0.75 (green solid curve) and fi' — —0.5 (blue dashed curve). The 
upper lines of the same colors represent the corresponding unextinguished source spectra. The best fit Reichart (2001) parameters for f3' = — 0.75(-0.5) are a source 
extinction Ay = 0.032(0.13), cj = 1.35(1.34), C3 = 8.1(2.7) and C4 = 30(6.9). The dip in the curves, between r* and g* , is the C3 redshfited 2200A bump. The steep 
drop through the u* observation is the C4 far-UV extinction («* samples the non-redshifted GRB spectrum from 1200Ato 1600A). The upper and lower dotted red 
lines correspond to an unextinguished f3' = —0.75 spectrum and the same spectrum extinguished by typical SMC-like extinction with Ay = 0.10. 




